Magnetic fluctuations in 2D metals close to the Stoner instability 



B. N. Narozhny and I. L. Aleiner 
Department of Physics and Astronomy, SUNY at Stony Brook, Stony Brook, NY 11794 

A.I. Larkin 

Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 55455 
L.D. Landau Institute for Theoretical Physics, 117940 Moscow, Russia 

We consider the effect of potential disorder on magnetic properties of a two-dimensional metal- 
lic system (with conductance 5 > 1) when interaction in the triplet channel is so strong that the 
system is close to the threshold of the Stoner instability. We show, that under these conditions 
there is an exponentially small probability for the system to form local spin droplets which are local 
regions with non zero spin density. Using a non-local version of the optimal fluctuation method we 
find analytically the probability distribution and the typical spin of a local spin droplet (LSD). In 
particular, we show that both the probability to form a LSD and its typical spin are independent of 
the size of the droplet (within the exponential accuracy) . The LSDs manifest themselves in temper- 
ature dependence of observable quantities. We show, that below certain cross-over temperature the 
paramagnetic susceptibility acquires the Curie-like temperature dependence, while the dephasing 
time (extracted from magneto-resistance measurements) saturates. 
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I. INTRODUCTION 

The itinerant ferromagnetism in metals has been the 
subject of extensive theoretical and experimental stud- 
ies. In clean systems the ferromagnetic instability is de- 
scribed by the Stoner criteriorocl, which defines the crit- 
ical value of the spin-exchange interaction constant (at 
which the system becomes unstable with respect to fer- 
romagnetic ordering). The value of the constant depends 
on the material. In most simple metals it is small (so that 
the electronic system exhibits paramagnetic response), 
while in palladium it approaches the critical valuea. The- 
oretically, the Stoner criterion is most easily obtained 
within the framework of-,the Hubbard modelnl or in Lan- 
dau Fermi liquid theoryBu. 

In two dimensions disorder tends to localize the elec- 
tronic systemEM However, if either the sample size L 
or the dephasing length L v is smaller than the localiza- 
tion length, the sample can still be considered as metal- 
lic. However, the interaction constants are renormalized 
from their clean values. In particular, the spin-exchange 
interaction constant was shown to flow towards strong 
coupling!! It is possible therefore to have a sample with 
the renormalized constant approaching the critical value 
while at the same time far from the Anderson localiza- 
tion. In this case a mean-field structure of the ferromag- 
netically ordered ground state was conjectured recently 
in Ref. H 

An important issue in the physics of disordered sys- 
tems is the role of mesoscopic fluctuations. Indeed, the 
Stoner criterion in its usual form reflects the tendency 
of the system to acquire uniform non zero magnetization 
(if the interaction constant happens to reach the critical 
value). However, for each realization of disorder the spin 
exchange interaction is non-local and random. The effec- 



tive interaction constant entering the Stoner criterion is 
essentially the interaction kernel averaged over the whole 
system. At the same time, averaging over some small 
part of the system might produce a value of the effective 
interaction constant different from the system- wide aver- 
age. In particular, some rare impurity configurations can 
lead to the locally averaged interaction constant in some 
region to satisfy the Stoner criterion, while the system- 
wide average does not. This would mean the appearance 
of non zero spin polarization in such regions. The similar 
effect in finite-size mesoscopic systems was considered in 
Ref. 0. 

In this paper we investigate the plausibility of such 
a scenario and its effect on magnetic and transport 
properties of the system. We consider a good metal 
(characterized by a large dimensionless conductance 
g = 2-Kh/e 2 Ru 3> 1) in 2D close to the instability, but 
still in the paramagnetic phase (at the mean field level), 
so that the naive mean-field value of the total magneti- 
zation of the system is zero. We describe the spin ex- 
change interaction by the effective averaged constant Fq 
with all the renormalizations already included (as shown 
in Ref. |^) and by the random, sample specific non-local 
susceptibility, which, when averaged over the area of a 
small region, gives the local effective interaction constant. 
Using the "optimal fluctuation" methodO we find that 
there is exponentially small but non zero probability to 
find such region with non zero spin polarization which 
we call a local spin droplet (LSD). With the exponen- 
tial accuracy this probability does not depend on the 
size of a LSD, therefore droplets of all sizes (up to the 
size of the order of the thermal length Lt — y/hD/T, 
(D is the diffusion constant and T is the temperature) 
can appear. The total spin of a LSD is also indepen- 
dent of its size. The effective interaction between LSDs 
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is determined by the correlations of the same non-local 
susceptibility. Since LSDs are extended objects, the cor- 
relation functions which determine the effective interac- 
tion constant have to be averaged over the area of both 
interacting LSDs. The oscillating parts of the correlation 
functions (which usually lead to the RKKYt^l interaction) 
do not survive this averaging. Instead, we find that the 
average value of the effective interaction is ferromagnetic 
and decays with the distance between LSDs only as a 
power law. However, at large distances fluctuations of 
the effective interaction constant exceed the average and 
the sign of the interaction becomes random. 

The contribution of LSDs to physical observables man- 
ifests itself at low temperatures. Indeed, as any sys- 
tem of weakly interacting moments at temperatures, 
higher than the point of magnetic ordering, the system 
of LSDs exhibit the Curie-like susceptibility. We show, 
that at not so low temperatures, this contribution ex- 
ceeds the (temperature-independent) Pauli susceptibility 
of the electron system. Likewise, the LSDs contribute 
to the dephasing time t v (extracted from magnetoresis=j 
tance measurements, see Ref. 0). Only, while usuallyUrtj 
the dephasing time in two dimensional electron systems 
behaves like t" 1 ~ T, the contribution of LSDs is tem- 
perature independent. Thus at temperatures lower than 
certain cross-over temperature the dephasing time satu- 
rates. The cross-over temperature is roughly the same 
for both quantities. When T — * 0, interactions of LSDs 
with each other or with itinerant electrons should lead 
either to the screening of the local spins or to forming of 
some spin glass state (due to the randomness of the in- 
teraction). We do not consider such regime in this paper. 

The paper is organized as follows. In Section || we re- 
view the basic physical description of the Stoner instabil- 
ity, establish notations and, in subsection [II C| we outline 
the effect of mesoscopic fluctuations. Next, in Section [II 
we disc uss th e formation of LSDs qualitatively. The sub- 
section [V A is devoted to the calculation of the probabil- 



ity to form a LSD. Then, in subsection IV B we find the 



distribution of the total spin oi-the LSD by means of the 
"optimal fluctuation" methodcil, which we adapt to the 
non-linear problem. Section |v| describes the interaction 
of LSDs, and in Section [Vl] we discuss the contribution 
of LSDs to physical observables. Our results are summa- 
rized in Conclusions. Some mathematical details of the 
non-local susceptibility correlations are relegated to the 
Appendix. 



II. STONER INSTABILITY 

The purpose of this Section is to recall t he b asic idea s 
leading to the Stoner criterion (subsections II A and [II B| ) 
and to contrast the situation in cle an sy stems to that in 
disordered metals. In subsection |II C| wc demonstrate 
the role of mesoscopic fluctuations and discuss how they 
can lead to formation of LSDs. For further details on 



Stoner ferromagnetism the reader is referred to the stan- 
dard textbooks Ref. H|| 



A. Renormalized paramagnetic susceptibility 

The paramagnetic response of a system of non- 
interacting electrons is described by the Pauli suscep- 
tibility which depends only on the electronic density of 
states at the Fermi level v. This can be seen from the 
following observation. The electron energy e enters all 
thermodynamic functions in combination e — \i with the 
chemical potential /i. The interaction energy of the elec- 
tron spin s with the external magnetic field h (which is 
proportional to —sh) can thus be considered as a shift of 
the chemical potential. Since it is proportional to s, the 
number of electrons which spin is aligned with h exceeds 
the number of electrons with the opposite spin, resulting 
in the total magnetization proportional to vh. 

Although quite general, the above argument relies 
on the fact that weak magnetic field does not change 
the energy spectrum of the electron system. Taking 
into account the electron-electron interactions, however, 
changes the distribution function of electrons and thus 
the electronic energy spectrum. As a result, some phys- 
ical quantities become renormalized from their bare val- 
ues. In particular, the paramagnetic susceptibility is 
renormalized by the the spin exchange interaction!] (we 
choose the units with Bohr magneton equal to unity) 
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1 + Fn 



(1) 



Here the parameter Fo is the effective dimensionless cou- 
pling constant of magnetic interaction between electron 
spins. Within, the phenomenology of the Landau Fermi- 
liquid theoryO it can be obtained by averaging the spin- 
exchange part of the Landau function over the Fermi sur- 
face. In the case when Fo < 0, the interaction tends to 
align the electrons spins, competing with the Pauli ex- 
clusion principle. If the interaction is strong enough, the 
gain in the magnetic energy exceeds the kinetic energy 
cost needed to realign the spins and the ground state of 
the system changes to the one with non zero total spin 
- it becomes ferromagnetic. According to Eq. (|l|) the in- 
stability occurs n wJien Fq = —I, which is known as the 
Stoner criterionBn 

More formally, the full susceptibility tensor is given 
by the commutator of the spin density operators a a (x, t) 
(hereafter x% denotes the two-component coordinate vec- 
tor) 



Xa/3(xi,x 2 ;ti - t 2 ) 



= id(ti - h)([a a (xi,ti),a0(x 2 ,t2)]), 



(2) 



where a, /3 = x, y, z. In the paramagnetic state of an 
isotropic system the tensor \ap is diagonal and isotropic 
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and can be expressed in terms of the transverse suscep- 
tibility x, 

which is determined in terms of the commutator of the 
spin raising and lowering operators a + and a~ similarly 
to Eq. (§). 

The transverse susceptibility can be evaluated in 
the generalized Hartree-Fock approximatioru, which 
amounts to summation of ladder diagrams in the particle- 
hole channel. In the Galilean invariant system the sus- 
ceptibility depends only on the coordinate difference and 
in the momentum representation is given byo 



n(g,w) 



l + UU(q,cj) 



(3) 



Here Tl(q,ui) is the electron polarization operator, which 
represents the susceptibility of the non-interacting elec- 
tron gas. In the limit uj — and q = it gives the 
Pauli susceptibility [since H(q — 0,u — 0) = v\. The 
parameter U is the spin exchange coupling constant. In 
the context of the Hubbard model it appears as a phe- 
nomenological parameter of the Hamiltonian. In the mi- 
croscopic Fermi-liquid theory it is the spin exchange part 
of the vertex function Y w averaged over the Fermi surface 
(compare with the scattering amplitude T2 in Ref. |^). 
The instability in the ground state of the system corre- 
sponds to the singularity in the static limit of the re- 
sponse function x{li^ — 0). The instability criterion is 
thus UH(q,u! = 0) = —1. At q = this corresponds to a 
tendency of the system to acquire spontaneously a uni- 
form (or ferromagnetic) spin density. The criterion for 
this instability 



Uv 



-1. 



(4) 



corresponds to the Stoner criterion if one identifies the 
Landau parameter Fq with the spin exchange coupling 
constant. 



be found by minimizing the thermodynamic potential Q. 
In the limit q = (uniform spin density) the minimum 
condition is 



(1 + F + Ba 2 )d = 0. 



(6) 



Here we have neglected the gradient term relative to the 
linear term, since the coefncientE-3 a ~ kp 2 and under 
our assumptions, g ^S> 1 and <?(1 + Fq) 1 (the latter 
assumption will be elaborated on in Section HI). There- 



fore the linear term dominates, 1 + Fq > kr/k 2 F even 
for large momenta k ~ l/l (I is the mean free path). 
As usual, in order to have a non-trivial solution a =/= 0, 
one should have 1 + Fq < 0, which is again the Stoner 
criterion. In the ferromagnetic phase the solution to 
Eq. (^) gives the total value of the induced spin density 
a 2 = -(l+F )/B. 



C. Mesoscopic fluctuations 

In the presence of disorder the local spin density a(x) 
depends on the particular impurity distribution and be- 
fore averaging can be taken as random. In the paramag- 
netic susceptibility (^|) the coupling U is determined by 
small distances and thus does not depend on disorder. 
On the contrary, the polarization operator Tl(xi, x 2 ,u>) 
includes large distances (since we are interested in the 
limit q = 0) and is thus strongly affected by disorder. 
Consequently, the paramagnetic susceptibility is random. 
Moreover, it depends on both coordinates since transla- 
tional invariance is lost and it is also non-local. However, 
if disorder is not too strong, then the non-local, random 
part of the susceptibility can be separated from the uni- 
form term, which is independent of disorder and repre- 
sents the susceptibility of the clean system. Preserving 
the form of Eq. (Q), the susceptibility of the disordered 
system can be found as a solution to the integral equation 



d 2 x 2 [(l + F )S(x 1 - x 2 )+F 1 (x 1 , x 2 )]x(x2, x 3 ) 



B. Magnetization energy 

The transition to the ferromagnetic state can also 
be described similarly to the Landau description of the 
phase transitions with the induced spin density as the 
order parameter. Close to the instability point, the ther- 
modynamic potential Q of the system can be expanded 
in powers of the spin density 3 




(5) 

where B > 0. The value of a of any particular state of the 
system (as described by the constants Fq, a, and B) can 



= v5(x 2 - x 3 ), (7) 

where Fi(xi,x 2 ) is a random quantity with zero mean. 
Similarly, the thermodynamic potential (||) becomes 

ft = J d 2 x(l + F )a 2 (x) 

+-!- J d 2 x 1 d 2 x 2 Fi{xi,x 2 )[a{xi)(r{x 2 )^ (8) 
1 f 4 

+— / JJcPxj-B [{xj}] ^a(x 1 )a(x 2 )j^(x 3 )a(x i )y 
v •* i=i 

Again, we are interested in the limit q — > (neglecting 
the gradient term; see the previous subsection and an 
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estimate below). In the second order term we have ne- 
glected the possibility of the presence of spin-orbit cou- 
pling. While the spin-orbit interaction can be taken into 
account, its presence docs not affect our main results (see 
Section ^ for discussion) . In addition to the fourth-order 
term written in Eq. (|9j) there is a term with different 
spin structure, namely [3{x\) x 17(22)] [^(23) x (7(24)] 
(see Appendix for details). In what follows we assume 
the simplest spin structure for a LSD a = (0,0, a). In 
this case the additional cross product term vanishes and 
we can treat the spin density as a scalar. 
Similarly to Eq. (§) 

the minimum of can be found 
from the (now non-local) integral equation 



{l+F )a( Xl 



d 2 x 2 Fi(xi, x 2 )a(x 2 ) 



f 4 

/ Y\_d 2 XiB [{xj}} a(x 2 )a(x 3 )(T(x4 : ) = 0. 



(9) 



The coefficient B \{xj}\ in the thermodynamic potential 
Eq. ([}]) is also random. In this paper we take both F 1 
and B to be Gaussian random matrices, with the distri- 
bution, which in compactified notation is given by 



; [Fx , B] oc exp 



- / (Ft; B)K 



- 1 / Fx 
B 



(10) 



The Gaussian approximation is valid while the expres- 
sion in the exponent does not exceed the dimensipnless 
conductance g, where the log-normal tail appearstJ. The 
integration is over all the variables of F% and B. The dis- 
tribution w[Fi,B] should be understood in the operator 
sense and will be used to evaluate the functional inte- 
grals below. The weight operator K. is constructed from 
the correlators 



K _ f (F1F1) 
' (BFi) 



(FiB) 
(BB) 



(11) 



which are discussed in detail in Appendix, where we give 
their explicit form. 

For a given realization of disorder, the equation (|^) 
might allow for some non-trivial solution a^(x). To es- 
timate the value of the total spin corresponding to such 
a solution, we write the spin density as 



r(°) 



(x) = (TQlpix), 



(12) 



where ip(x) is normalized to unity, therefore both ^(x) 
and (To have dimension of inverse length. The total value 
of the spin is determined by 



S = <7q / d 2 xip(x), 



(13) 



while (To can be found from Eq. (^|) (in the case when it 
allows a non-trivial solution) in the form 



1 + F + F{ 



(0) 



(14) 



Here the constants F^ and B^ are the "matrix ele- 
ments" of the non-local operators Fi(xi,Xa) and B [{xj}] 
(if ip(x) is interpreted as a "wave function" ) 



?(o) _ 



J d 2 xid 2 x 2 F 1 (xi, x 2 )ifj(xi)ip(x2) (15a) 



4 4 
f\{d 2 x l B[{x ] }]\{^{x k ). (15b) 

^ i=l k=l 



Consider now a metal close to the Stoner instability, 
so that < 1 + F <C 1. While, as follows from Eq. (0), 
the averaged, uniform spin density is zero, solutions ( |14| ) 
to the non-local equation (|9|) might exhibit non zero spin 
density in some rare regions, where due to a particular 
configuration of impurities the non-local part of the sus- 
ceptibility F\ is negative. If there are several regions with 
non zero total spin, then the spin-spin or magnetic inter- 
action between them will contribute to the ground state 
energy of the system and correspondingly to the mag- 
netic susceptibility. If such interaction favors some kind 
of ordering of the spins, then the appearance of these re- 
gions can change the magnetic response of the system, in 
other words change the ground state. 

Fluctuation effects in systems close to a phase tran- 
sition have been studied extensively (see, for instance, 
Rcf. [l4|,[l(] |l8|). In particular, the picture of smearing 
the transition point by formation of fluctuation regions 
with non-zero value of an order parameter was consid- 
ered in Ref. |1| In order to determine a value of the 
order parameter in a fluctuation region a solution of the 
non-linear Ginsburg-Landau equation in the presence of 
disorder was needed. Two issues make this case different 
from ours. First, in Ref. |l^ the fourth-order term B was 
not random. Second, the fluctuations of the order pa- 
rameter were assumed to be local. Therefore, the white 
noise approximation (i.e. approximating the correlatino 
functions (|ll]) by the delta-functions) for the disorder 
was appropriate. As a result, the fluctuation regions did 
not interact and the percolation scenario of the phase 
transition was needed. In our case the non-locality of 
fluctuations (expressed in terms of the correlation func- 
tions ( p"l| ) leads to interaction between LSDs which re- 
sults in a change of behavior of observable quantities as 
discussed in Section VI. 



III. QUALITATIVE DISCUSSION 



In the previous subsection |II C| we indicated how LSDs 
- local regions with non zero spin polarization - could 
appear in a metal close to the Stoner instability due to 
fluctuations in impurity distribution. Here we estimate 
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qualitatively the probability to find a LSD and the value 
of its total spin. 

Treating a LSD as an open region of the size R we 
can characterize it by the Thouless energy Et = D/R 2 
(where D is the diffusion constant). The inverse of the 
Thouless energy is the "escape time" r esc — E^ 1 , which 
is the time it takes for the diffusing particle to leave the 
LSD. This time scale serves as the infrared cut-off for 
the correlation function, which describes tlie mesoscopic 
fluctuations of the density of states (DoS)t£l 



(p(e)p(e + ^))-Re 



En 



R 2 d 2 Q 



-Re- 



{—iuj + DQ 2 
1 



"^2 



-iuj + T e . 



(16) 



The magnetic energy of the LSD, written in terms of its 
total spin S is 



E(S) = 6 1 (l + F )S 2 +6E(S), 



(17) 



where Si is the mean level spacing and SE(S) denotes 
contribution of all non-linear terms in Eq. (0). All the 
terms contributing to SE(S) are random and can be ex- 
pressed in terms of the fluctuating (random) DoS 



<5iS s i 

SE(S)= J dsx J ds 2 [p{s 2 ) + p(-aa)] 



(18) 



The averaged SE(S) equal to zero, but the average of 
its square ((8E(S)) 2 ) is not and it is determined by the 
correlator ( jig ) 



((SE(S)f 



a 



5> 5 . 



(19) 



In this paper we restrict ourselves to consideration of 
metals close to the instability point, where the overall 
spin of the LSD is small S <C g. In this case we can treat 
6E(S) as a Gaussian random quantity. Moreover, we can 
expand it in powers of S/ g so that the magnetic energy 
of the LSD becomes 



E(S) =5 1 (l + F )S 2 + ( 
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g 3 



(20) 



where £ is a random Gaussian variable with the distribu- 
tion 



P(C) ex cxp (-C 2 ) • 



(21) 



Note that there is only one random quantity ( in Eq. ( |2l[ ) . 
The equation (|o|) is valid only for S <C g. That is why 
energy minima at£ > are spurious and should not be 
considered. 

To find the distribution of the spin value, we need to 
minimize the energy E(S). Differentiating Eq. ( p0| ) we 
obtain the equation 



A S" 3 
1 + F + ± )S - C^r = 0. 

gj g 6 



(22) 



This equation allows for non-trivial solutions S > when 
C < — g(l + Fq). As was noted above, we are working in 
the regime where g(l + Fq) ^> 1, so that the fluctuation 
that creates the LSD is rare indeed. To determine the 
distribution of spin values, we solve the equation (|22j) for 
C and substitute in the Gaussian distribution Eq. (|2l|). 
As a result, we estimate the distribution (up to numerical 
coefficients) 



P(S) = cxp - 



(i 



(23) 



This distribution is only valid when g(l + Fq) 3> 1 and 
S <^ g, therefore we can expand the denominator in the 
exponent without loss of accuracy 



P(S) = exp [-.g 2 (i + F ) 2 - (1 + F ) 2 S 2 ] 



(24) 



Integrating over the spin S we estimate the probability 
to find a LSD 



Fcxexp^g^l + ^o) 2 ) 



(25) 



It is determined by the first term in the exponent 
Eq. (|24|). The second term determines the typical value 
of the spin of the LSD 



(26) 



l + Fn 



Remarkably, this value and the probability Eq. (|25|) do 
not depend on the size R of the LSD, which is the main 
qualitative result of this Section. Note ihat it is similar 
to the result for zero-dimensional grainst! Also the spin 
value is independent of the dimensionless conductance g. 
These facts determine the contribution of LSDs to the 



physical observables considered below in Section VI 



The applicability of the consideration of this Section 
and quantitative results of Section [V is limited by two 
requirements. First, the probability Eq. (|2^) must be ex- 
ponentially small, so that g(l + Fq) 1. Second, the 
Gaussian approximation Eq. ( ^l| ) [similarly to the dis- 
tribution Eq. ([h])] is valid while the expression in the 
exponent is smaller than the dimensionless conductance 
of the systemta, so that in Eq. (||) g 2 (l + F ) 2 < g. 
Combining the two limits, we obtain the region of appli- 
cability of the results Eqs. (23) - (|26|) as 



l/ 5 « 1 + Fq < l/y/g. 



(27) 
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IV. LOCAL SPIN DROPLETS 

In this section we calculate the probability to find a lo- 
cal region with non zero spin, which we call a local spin 
droplet (LSD) and the value of the total spin of the LSD. 
The spin and the spatial profile of the LSD can be found 
from the non- linear equation (m . In subsection IV A 



we 



show that to the exponential accuracy the probability to 
find the LSD Eq. ( p5| ) is captured by the linear part of 
Eq. (^|), while the no n-line ar term fixes the spin value, as 
shown in subsection IV B. 



A. Probability to form a LSD 

The calculation of the probability to find a rare fluc- 
tuation leading to formation of a LSD can be performed 
along the lines of the argument used in Ref. [TJ to calcu- 
late the exponentially small tail in the density of states 
(DoS) of a particle in a random potentially. In the quan- 
tum mechanical problem, considered in Ref. O, one looks 
for such fluctuation of the random potential that creates 
a low energy bound state, thus leading to non zero DoS at 
that energy. The probability to form the bound state is 
determined by the distribution of the matrix elements of 
the random potential. While being exponentially small, 
the probability should be maximized by choosing the "op- 
timal" fluctuation of the potential. 

To gain some intuition about how the optimal fluctu- 
ation method can be applied to the problem at hand, in 
this section we consider the linear part of the equation 
(^), disregarding for a moment the higher order B term. 
Such an approach can be justified by observing that close 
to the instability the non-linear term in Eq. (^|) is small 
compared to the linear ones, since the induced spin den- 
sity on average is equal to zero. The non-linear term 
stabilizes a non-trivial solution and fixes its amplitude, 
while the existence of such a solution can be uncovered at 
the level of the linear problem. Thus the linear equation 
captures the main contribution to the probability and at 
the same time demonstrates the similarity of our problem 
to the problem of tails in the DoS as well as the peculiar 
differences. 

We can write the linear equation in the operator form 



(28) 



We write ip(x) instead of a(x) to stress the point that 
the linear equation does not allow us to determine the 
value of the spin but only the spatial profile of the LSD. 
Therefore, Eq. ((2^) is simply the eigenproblem for the 
operator F\ and as such does not fix the normalization 
of eigenfunctions ip(x), which we are free to normalize 
to unity for convenience. Since the eigenvalue problem 
Eq. ( p8|) i s similar to the quantum mechanical problem 
of Ref. |11|, we can adopt the language of the Schrodinger 
equation, with the (now integral) operator F\ playing the 



role of the "Hamiltonian" , E the "energy" and ip(x) the 
wave- function. 

For some particular realizations of the random poten- 
tial in Eq. (|28|) there is a low energy bound state with 



the energy iJ [{-Fi}] = F{ ' [given by Eq. ( |15a| )l, result- 
ing in non zero DoS at this energy. For energies close 
to £?o[{-Fi}] only the bound state contributes to the DoS, 
which before averaging over disorder is given by the single 
delta-function 



p{E) = 5{E - E [{F X }]). 



(29) 



Averaging this DoS over disorder takes into account con- 
tributions of all possible realizations of the random po- 
tential leading to such bound states and thus results in 
an exponentially small but smooth function of energy E. 
This function is proportional to the probability to find 
the bound state at energy E. In particular, for the spe- 
cial value E = — (1 + Fq), it would give the probability to 
find the non-trivial solution to the linear part of Eq. (|^) 
(or to find the LSD). 

The random quantity in the linear problem Eq. (p8|) is 
Fi itself. Its distribution is obtained from Eq. (Jll|). The 
averaged probability is then 



(p(E)) = J V[F X ]5{E - E [{F X }]) 



x exp 



d X\d x 2 d y\d y 2 A 



(30a) 



A = F 1 (x 1 ,x 2 )K FF [{x j },{y j }}F 1 (y 1 ,y 2 ), (30b) 

where K FF [{xj}, {yj}] is the inverse of the correlator 
Eq. (]A|) i.e. 



d 2 y 1 d 2 y 2 K FF [{ Xj }, {%}] K FF [{%}, {z 3 }} 



5{x\ - zi)8{x 2 - z 2 ). 

In the "optimal fluctuation" approach one has to eval- 
uate the integral Eq. (^) in the saddle point approxima- 
tion. To find the saddle point one has to minimize the 
exponent A of the Gaussian probability with respect to 
all functions Fi(x±,x 2 ) subjected to the condition 



J B = F [{F 1 }] 



(31) 



represented by the delta- function in Eq. (|3C|). This in- 
volves solving the equations 

[A + XE [{F 1 }]] = 0, 



5Fi(xi,x 2 [ 



where A is the Lagrange multiplier to be found from the 
condition Eq. (|3l|). 

The saddle point solution for Fi represents the optimal 
fluctuation of the random potential, 



G 



Fi(x 1 ,x 2 ) 



A / d 2 z 1 d 2 z 2 K FF [{ay}, { Zj }] ip{zi)iP(z 2 ), (32) 

given in terms of the eigenfunction ip(x) corresponding 
to the eigenstate E . Substituting Eq. (|32|) into Eq. ( |2"§| ) 
we obtain 

Eip{xi) = A / d 2 x 2 d 2 zid 2 z 2 tp(x 2 ) 



xtf FF [{^},{^}] ^1)^*2), (33) 

which, together with the normalization condition 
J d 2 x\ip(x)\ 2 = 1, constitute the analogue of the non- 
linear Schrodinger equation of Ref. [ll]. 

Both the saddle point value of A and ip(x) should be 
found from the non-linear equation (^3|) and the nor- 
malization condition. The averaged DoS is given by the 
Gaussian probability in Eq. ([50|) , evaluated at the saddle 
point, 



(p(E)) ~ exp 



E 2 



2I FI 



(34) 



where 



iff 



d z\d z 2 d x\d x 2 



xip(xi)ip(x 2 )K F F [{xj}, {zj}] ip(z 1 )^(z 2 ). (35) 

The integral Iff is dimensionless (independent of any 
length scale) since we require the eigenfunction tp(x) to 
be normalized, so that ip(x) ~ where R is roughly 
the size of the LSD. As we are discussingthe single 
LSD, all four of the eigenfunctions in Eq. (J35|) are cen- 
tered around approximately the same point, therefore R 
is the only scale in Eq. (pq). In this case the depen - 
dence of the correlator Kff on R is given by Eq. ( A12 ), 
Kff oc g~ 2 R~ 4 . Thus the integral Eq. ( p5| ) and the 
probability Eq. (|34| ) are independent of the size of the 
LSD R. 



The probability to find the LSD is given by Eq. (34) 
evaluated at the point E = —(1 + F ), 



P ' 



,-7ff 2 (l+-F ) 2 



(36) 



where 7 is the numerical factor which is given by the di- 
mensionless counterpart of the integral Eq. (|3^). This 
result is valid while the number in the exponent is large, 
which corresponds to the lower limit of applicability of 
our consideration 1 + F ^> 1/g (see the last paragraph 
of Section H|). 

To determine the numerical coefficient 7 in Eq. (|3^) we 
need to know the precise form of the eigenfunction ip(x). 
We could not solve the non-linear integral equation ( |33| ) 
analytically and used a variational approach. Since the 



kernel Kp f in the integral Eq. (35) decays as R~ A at 
large distances, the optimal ip( x ) is a limited-range func- 
tion (i.e. its normalization integral is determined on a 
limited interval of x). To estimate the upper limit of 
7 we take the variational function of the Gaussian form 
xj;{x) = n- 1 / 2 R- 1 exp(-x 2 /R 2 ), substitute in Eq. (||) 
and evaluate the integral numerically. The resulting es- 
timate is 



1 



FF 



0.8ttV 



(37) 



and the numerical factor in Eq. ( |3q ) is thus 7 ~ 3.9. 

The most important feature of the result Eq. (|3^) is 
its independence on the size R of the LSD in agreement 
with the qualitative results of Section III. This means 
that LSDs of small and large sizes can appear with equal 
(in the exponential sense) probability, given the suitable 
fluctuation of the impurity configuration. We shall return 
to this point below, when we consider the interaction be- 
tween LSDs. 



B. Spin distribution 

The argument that led us to the probability to find an 
LSD Eq. ( |36"|) is not complete because it does not help us 
to determine the value of the spin of the LSD. This fol- 
lows from our consideration of the linear part of Eq. (j^) 
only. To determine the value of the spin we must solve 
the full non-linear equation. 

The single-mode approximation to the solut ion o f the 
non-linear problem was outlined in subsection II C . The 



formal solution for the amplitude of the spin of a LSD 
is given by Eq. (|l4|). Again, as we did above in the case 
of the linear problem, we employ the optimal fluctuation 
method to find the probability to form a LSD, which is 
characterized by the total spin S. However, this time our 
task is simplified since the probability to form the LSD 
(regardless of its spin) has already been found. We now 
need to find how that probability depends on the spin 
value S. Therefore we take the function ip(x), which 
describes the spatial profile of the LSD from the linear 
problem and focus on the distribution of the spin ampli- 
tude a. 

Similar to Eq. (p9|), the probability to find an LSD 
characterized by the spin density amplitude a is given 
by the delta-function 



p(a 2 )=S(cx 2 -a 2 [F 1 ,B}), 



(38) 



which we write in terms of a 2 for convenience. 

Averaging over disorder is performed as it was done 
for the linear problem [see Eq. (|30|)]. Only now we have 
two random quantities, F\ and B, therefore we need to 
average with the distribution Eq. ( |l0| ) 



(p{° 2 )) 



V[F 1 ,B]8{a 2 - al)w[F u B], 



(39) 
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To find the saddle point (or the optimal fluctuation) we 
have to minimize the exponent 



Au = A <7 



2 , l+Fo+F± 



(0)' 



(F 1 B)^- 1 



Fx 
D 



(40) 



where F^ and £?(°) are the integrals Eq. JTe| ) and A is 
again the Lagrange multiplier. The saddle point equa- 
tions are given by 



Kr 



i (Ft 
B 



A 



, (41) 



where the functions ip( x ) appeared after differentiating 
the integrals Eq. (|l5|) with respect to F% and B. Mul- 
tiplying both the left-hand and the right-hand sides of 
Eq. (|d|) by K. we obtain two equations with F\ and B 
in the left-hand side only. We then multiply the first 
equation by two functions ip(x) and the second by four 
functions ip( x ) an( i integrate over their variables to ob- 
tain the algebraic equations, with integrals Eq. ( |l5| ) as 
the unknowns 



F^ = 



2B(°) 



IFF 



+ o»I 



FB 1 



* (0) = ^W(^+^) 



(42a) 



(42b) 



where Iff is given by Eq. ( |35| ) and Ifb and Ibb are 
similarly defined as 

Ifb = J d 2 xid 2 x 2 d 2 y i ip(x 1 )ip(x 2 ) 

i—l 

4 

xKFB[{xj},{yj}}T[i>(yk), (43a) 



fc=l 



Ibb = / JJd 2 a;jd 2 yi V'fafc) 
i=i fc=i 

4 

x^bbK^}, {»,-}] JJ^(W)- (43b) 



fc=i 



Equations (f42| ) can be easily solved and we have for the 
optimal fluctuation 



f (0) = Vx- Iff + ct2/fb 



V2 {Ifb + (J 2 IbbY 
B (0) = iVA^ OfbT^JbbT. 



(44a) 



(44b) 



Now we only need to find A from the constraint a 2 = 

-(1 + F Q + F^ 0) )/BW or, equivalently, to find the sad- 
dle point solution for A. Substituting the saddle point 
solutions ( p4| ) into the constraint, we find 



s/\ = -(l + F ) 



FB 



IBB 



Iff + 2a 2 I F B + <J 4 Ibb ' 



(45) 



The exponent A n i at the saddle point Eq. 
by 



J) is given 



A n i (A) = 



A 



2£(°) 



a 2 S(°) + P[ 0) + 2(1 + F ) 



(46) 



Finally, evaluating the exponent ( flq ) for the optimal 
value of A Eq. (f45|), we obtain the probability to find 
an LSD with the spin density amplitude a 



(p(a 2 )) ~ exp - 



(l + ^o) 2 



2(I FF + 2<7 2 Ifb + v 4 Ibb) 



(47) 



which for a — coincides with the result Eq. (^) of the 
linear problem, just as the result Eq. (|^) of the qualita- 
tive argument above. 

The final step is to rewrite the distribution Eq. (|47j ) 
in terms of the spin S of the LSD. Converting the spin 
density amplitude a in Eq. (^7]) to the spin S by means 
of Eq. (|l3|), we obtain the final expression for the spin 
distribution 



p(S 2 ) ~ exp 



7 ff 2 (l + Fo) 2 
l-2af£ +/?{ 



where the numerical factors are 
9 2 \Ifb\ 

a = 7y 

Iff [fd 2 x^(x)Y 
a 9 4i bb 



Iff [j (Pxx/j(x)]' 



1.92, 



7.36, 



1 



2g 2 I FF 



3.9. 



(48) 



(49a) 



(49b) 



(49c) 



The factor 7 is the same as in the linear problem and is 
listed here for completeness. 

The factors Eq. ([l9|) are now independent of the size 
R of the LSD. To see that, one needs to notice that all 
the correlators fC depend on this scale in the same way 
(when a single LSD is considered, so that R is the only 
scale in the integral) JC ~ i?~ 4 [see Eq. ( A16Q ]. The wave 
function tJj(x) is inverse in R. Therefore, the integrals 
Eq. (43), unlike the integral Eq. (|35|), which appears in 
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the solution to the linear problem, do depend on R, since 
they contain different number of functions ip{x) 



R 2 

Ifb oc —, 

9 



IBB OC —~. 

g 6 



(50a) 



(50b) 



The different size dependence of the integrals Eq. ( |35| ) 
and I FF is compensated in Eq. (49) by additional fac- 
tors of J d 2 xip(x). 



The distribution Eq. fl48| ) is thus the same as the quali- 
tative result Eq. (23) only now with coefficients Eq. (49). 
The coefficients were evaluated numerically using the so- 
lution ip(x), which follows from the linear problem. The 
coefficients are positive [the negative sign of Ifb being 
taken into account explicitly in Eq. (Eq)]. The coeffi- 
cient /3 is positive, ensuring convergence of the expansion 
Eq. (|2^), which again is the way to interpret the distri- 
bution Eq. (jig 1 ). In the calculation leading to Eq. jjq 1 ) 
the limitation to small spins follows from the separation 
procedure Eq. (|l2| ) since it is valid only close to the in- 
stability, where typical spins are small. 

The typical value of spin of the LSD is still given by 
Eq. (ptf). The spin S turns out to be large, 1 <C S <C g, 
so on length scales larger than R or at high enough tem- 
perature, LSDs behave as classical moments. However, 
LSDs are extended objects and can have any size with 
equal probability (with the exponential accuracy). Thus 
their spins can not be considered as local moments, es- 
pecially when discussing their interactions. 



V. INTERACTIONS BETWEEN LSDS 



The importance of LSDs is that their appearance can 
dramatically change the magnetic response of the system. 
At high enough temperatures, we can consider them as 
independent, classical moments, thus we expect the sys- 
tem to be the usual paramagnet with the susceptibility 
described by the Curie law (but with Curie constant dif- 
ferent form that of the free electron gas). As the tem- 
perature becomes smaller, the system might change its 
ground state in a way that depends on the interaction be- 
tween LSDs, in particular on its sign and typical range. 




R 

FIG. 1. Two LSDs located around points a and b and the 
three averages of the non-local susceptibility Fi. The two 
averages and F^ describe correlations within each sin- 

f 12) 

gle LSD, while F£ describes correlations between different 
LSDs and thus determines the interaction 



Let us first recall the basic physics of local moments in 
a metallic system. In a system of local moments there 
are two competing types of,interaction. First, there is the 
direct contact interactional, which sign does not change 
with distance between the moments, but the amplitude 
decays exponentially beyond the correlation length. This 
interaction tends to turn the system into a ferromagnet. 
Then there is the RKKY interactions, which as a func- 
tion of distance oscillates and decays only as a power law 
(~ R~ 2 in 2D). The RKKY interaction between local mo- 
ments tends to form a spin glass at low temperaturestD. 

Our case is different. As we have shown in this pa- 
per, in a metal close to Stoner instability there is a non 
zero probability for LSDs to be spontaneously formed. 
This probability is independent of the size of LSDs, so 
that droplets of all sizes can appear. Therefore interac- 
tion between the LSDs can not be described in the same 
way as interaction between local moments. Rather, it is 
given by the non-local susceptibility F\(xi,X2), averaged 
over the area of the interacting LSDs. Calculation of the 
average, previously denoted as , involves integration 
of F\{x\, X2) with two functions tp(x — a) which describe 
the spatial profile of the LSDs. In this section we use 
the same Gaussian functions we used to evaluate the in- 
tegrals Eq. ( |35| ) and Eq. ([43|). Only now these functions 
carry explicitly the dependence on the coordinate a of 
the center of the LSD. 

Consider now two LSDs separated by distance L much 
larger than the size of both LSDs L ^> R, see Fig. 1. 
In this argument we take both LSDs to be of the same 
size R, but it can be easily generalized for the case where 
the interacting LSDs differ is size substantially. Choos- 
ing different combinations of the wavefunctions ip(x — a) 
and ip(y — b) describing the two LSDs, we can form three 
different averages 

F^ — J d 2 x\d 2 X2Fi(x\, X2)'ip(xi — a)ip(x2 — a), (51a) 
F[ 2) = / d 2 y 1 d 2 y2F 1 (y 1 ,y2Wy 1 - b)^(y 2 - b), (51b) 
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(12) _ 



d 2 xd 2 yF 1 (x,y)il}(x - a)ip(y - b). 



(51c) 



These averages, as random quantities, have a Gaussian 
distribution with the weight determined by the correla- 
tion function (F\Fi). We now assume that the LSDs have 
been already formed. The saddle point solution F\ [see 
Eq. (p4)1 describes a single LSD, while the interactions 
are determined by small deviations from the saddle point. 
Therefore the weight of the distribution of the averages 
Eq. (|5l|) can be evaluated at the "optimal fluctuation" 
point Eq. (|32|). Then the distribution can be written as 



exp 



F 



(i) 



F 



(2) 



F (12) ) £-1 



\Fl 



(12) 



(52) 



Iff 


■h 






Ji 


Iff 


i). 


(53) 


■h 


J-2 







with the weight matrix 



C 



The elements of the weight matrix C arc obtained by av- 
eraging the correlation function Kpp over the area of the 
LSDs, i.e. integrating with four wavefunctions ip(x — a), 
similar to the integral Iff [see Eq. (35)]. The difference 
from the case of the single LSD is that now all but two 
of the elements of C depend on two different lengthes : 
the size of the LSD R and the distance L between them. 



Thus the estimate Eq. (A16) for the correlation function 
Kff does not apply. The elements Ji are given by the 
integrals 



■h = J d 2 x Y d 2 x 2 d 2 y\d 2 y 2 K F F [{%j}, {Vj}] 



xi/j(x 2 - a)ip(yi - b)tp(y 2 - b) 



0.15 R 4 
L 4 '' 



a) 



(54a) 



J-2 = J d 2 x 1 d 2 x 2 d 2 y 1 d 2 y 2 K F F [{Xj},{Vj}] 



xij){x\ - a)*l>(x 2 - b)ip(yi - b)ip(y 2 - b) 



The integrals Eq. (Q) were calculated in the leading or- 
der in R/L. Clearly, J 2 ^ J 3 exactly, but the difference 
comes in the numerical factor under the logarithm, which 
we here neglect. This does not have any bearing on our 
conclusions. The l/L 4 dependence of all the i ntegrals 
follows from the frequency integral in Eq. ( All ), which 
is determined by the Thouless energy corresponding to 
the largest length in the problem, which is now L. 

The integral Ip f is independent of all length scales and 
therefore is much larger than any of J,; . Thus the weight 
matrix £ -1 in Eq. (|52|) to the leading order in R/L is 



1 w O.SttV 




(55) 



where the dimensionless function f(R/L) is given by 



f(R/L) « 0.72£ In 2 R 



L 4 



L ' 



(56) 



To describe the interaction between LSDs we need to 



find the distribution of F n 



(12) 



under the condition that 



the two LSDs exist, namely that < and f[ 2) < 0. 
This is given by the conditional probability distribution 

W[F[ 12) ]= ± J V[F^,F^}6(-1 -F - fP) 



x6(-l -Fa- Fi 2) )w[F^',F^',F^'], (57a 



(!) i?( 2 ) p( 12 )l 



N 



V[F\ 



(Di 



-d(-l-F -F, 



(57b) 



where w is the distribution Eq. (|52]). Since the weight 
matrix in Eq. ( ^2| ) is a c-number, the integrals in Eq. ( p7] ) 
are usual Gaussian integrals and not functional inte- 
grals. The (9-functions in Eq. ( E7t ) make the exact in- 
tegration in terms of elementary functions impossible, 
but we can use the small parameter f(R/L) <C 1 to es- 

( 12) 

timate VFfi 7 ^ ] with exponential accuracy, which is all 
we need to describe interaction between LSDs. Up to the 
pre-exponential factor 



1 R\ 2 R 



hr 



6ttV L 4 L 



1 R A , 2 R 

In — . 



6tt 2 5 2 L 4 L 



(54b) 



J 3 = J d 2 xid 2 x 2 d 2 yid 2 y 2 K F F [{xj},{Vi}] 



xip(x 1 - a)ip(x 2 - b)ip(yi - a)ip{y 2 - b) 



(54c) 



W[F\ 



(12)l 



oc exp 



0.8ttV 



(f i (12) +2(1 + ^ )/ 



(58) 



so that the distribution is a Gaussian (as it should be 
since we considered F\ to be a Gaussian random quan- 
tity from the very beginning). 

The average f[ 12 ^ given by the distribution Eq. ( p8| ) is 
shifted from zero to the negative value —2f(R/L)(l+F ), 
i.e. the average interaction appears to be ferromagnetic. 
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However, the distribution Eq. ( |58| ) also allows for strong 

( 12) 

fluctuations of F± '. These fluctuations can be esti- 
mated as 



am 



(X2)x 2 \ 



1 



(12) \ 



(59) 



The function / [see Eq. (|56"|)] decreases with the distance 
between LSDs. Therefore at large enough distances the 

(12) 

fluctuations of F± exceed the average and the sign of 
F^ 12 ^ becomes random. 

The cross-over distance L* can be estimated as (from 
the condition {{5f[ 12) ) 2 ) / {f[ 12) ) 2 ~ 1) 



L* ~ i?v/s(i + ^b)> R. 



(60) 



In the model considered in Ref. |Ts| mesoscopic fluctua- 
tions become uncorrelated beyond the magnetic length. 
In our case the role similar to that of the magnetic field 
in Ref. ^ would have been played by the spin-orbit cou- 
pling, which we have so far neglected. However, taking 
the spin-orbit coupling into account does not change the 
main results of this paper. As we have shown above, the 
typical interaction between LSDs is random due to large 
fluctuations in Eq. ( p8[ ) . Should we include the spin-orbit 
coupling, we would need to compare the spin-orbit length 
to L* in order to determine the length beyond which the 
interaction becomes random. But since the typical dis- 
tance between LSDs is exponentially larger than L* the 
exact value of such cross-over length is not very impor- 
tant. 



This distance has to be compared with the typical dis- 
tance between LSDs. The latter can be estimated as fol- 
lows. The concentration of LSDs has to be proportional 
to the probability Eq. p6]) 



-79 2 (1+F ) 2 



(61) 



where R is the characteristic length of a LSD. There- 
fore the typical distance between LSDs is exponentially 
large in the parameter g(l+Fo) 3> 1, while the cross-over 
length L* is large only as a power of the same parameter. 
Thus, the interaction between typical LSDs has random 
sign. 

The energy of interaction between typical LSDs decays 
as the second power of the distance between them 



Ui 



1 



1 



- 7ff 2(l + F ) 2 



vgL 2 ugR? 



(62) 



where v is the density of states and we have estimated 
the typical distance between LSDs from Eq. (|6l|). 

The results of this Section resemble the results of 
Ref. [l7|, where the electron-mediated interaction between 
magnetic moments in two ferromagnets separated by a 
disordered metal was also found to have a random sign. 
The diference between our model and that of Ref. [I?] is 
that in our case the ferromagnetic regions (LSDs) are cre- 
ated by the very same impurity configurations as those 
responsible for the interaction between LSDs. Therefore, 
the interaction between LSDs can be considered only un- 
der the condition of their cxistance and thus is character- 



ized by the conditional probability distribution Eq. (58). 

(12) 

As a result, the average F± is negative, i.e. ferromag- 
netic, and dominates the fluctuations at distances smaller 
than L* . However, since the typical distance between 
LSDs is exponentially large, the typical interaction has 
random sign. 

This fact is not surprising, since we are dealing with a 
disordered system where one could expect to find some 
spin glass phase (at T = 0, similar to the case of a super- 
conductor in weak magnetic field, considered in Ref. [l^. 



VI. CONTRIBUTION TO PHYSICAL 
QUANTITIES 

In this section we estimate the contribution of the 
LSDs to observable physical quantities. We consider the 
examples of paramagnetic susceptibility and the dephas- 
ing time. To be observable, the contribution of LSDs 
should exceed the regular contribution of the electron 
system. We show that, for both quantities, it happens at 
the same temperature T*, below which the temperature 
dependence of both quantities changes. The Pauli sus- 
ceptibility crosses over to a Curie-type 1/T dependence, 
while the dephasing time saturates and becomes temper- 
ature independent. 

Since the typical interaction energy Eq. (|62| ) is expo- 
nentially small, LSDs behave as weakly interacting mo- 
ments and their contribution to the paramagnetic suscep- 
tibility is given by the Curie law xlsd = C/T, where the 
Curie constant C is proportional to the square of the spin 
S of the LSD [see Eq. (Pq )] and the density (or concen- 
tration) of LSDs Eq. (pThT We have calculated the prob- 
ability Eq. ( |36| ) with the exponential accuracy, therefore 
the pre-exponential factor in the density of LSDs should 
be given by the typical length characterizing the distribu- 
tion of LSDs. This dependence could be elucidated from 
the following dimensional argument. The probability to 
find an LSD in unit volume (or, rather, area, which we 
still denote by V) of size in between R and R + dR, with 
spin in between S and S + dS is given by the distribution 

dW = P {S 2 )dVdSd 2 R (63) 

Clearly, the dimensionality of p{S 2 ) is L~ 4 , where L is 
the characteristic length. Since the exponential in p(S 2 ) 
[see Eq. (48)] is independent of length, the size of a LSD 
can not be pinned to any length in the system. Since in- 
dividual LSDs contribute to Curie susceptibility indepen- 
dently, we need to sum over all possible sizes and given 
the L~ 4 dependence, the main contribution comes from 
LSDs of the smallest possible size, namely the mean-free 
path I. Thus, contribution to the susceptibility (up to a 
numerical coefficient) is 
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1 



XLSD OC 



TTg(l + F o y 



_ 7ff 2(l + F ) 2 



(64) 



where r is the scattering time. Here we used the typical 
value ( |I3| ) of the spin of a LSD. This contribution has to 
be compared with the Pauli susceptibility Eq. (0) , which 
is temperature independent. As a result, the contribu- 
tion Eq. (|^) dominates at temperatures smaller than T*, 
which up to numerical factors is 



physical quantities suggests that our scenario of the mag- 
netic fluctuations in a metal close to the Stoner instability 
can be experimentally observed. The LSDs lead to the 
saturation of the dephasing time at low temperatures. 
If such saturation is observed, one should look at the 
behavior of the paramagnetic susceptibility in the same 
temperature range. If the LSDs are present in the system, 
then the onset of the Curie-like temperature dependence 
should also be detected. 



T* 



rg(l + F ) 



-19 2 (l+F ) 2 



(65) 



This temperature has to be compared with the typical 
interaction energy Eq. (^). Substitutin g th e mean- free 
path I for the chracteristic size R in Eq. we find the 
ratio U t /T* to be 



(l + F ( 



oj 



T* ~ g(l + F ) 



< 1. 



(66) 



Therefore, the Curie behavior Eq. ( |64[ ) persists over a 
wide temperature range U t < T < T* . 

Similarly, we estimate the contribution of LSDs to the 
dephasing time t v . LSDs are regions where the impu- 
rity configuration makes it energetically favorable for the 
electrons to align their spins. Some other electron enter- 
ing such region will "feel" the overall polarization as if it 
were magnetic field. The corresponding dephasing time 
can be estimated as the square of the Zeeman energy di- 
vided by the Thouless energy (at the size of the LSD). 
More formally, since the interaction with the polarization 
stems from the exchange-interaction, we can find t v from 
the perturbation theoryfj 



1 



F 2 



ri 



V 



r 2 



V 



CT(r 1 )cr(r 2 ))X>(r 1 - r 2 ) 



where T>(r\ — r 2 ) is the (electron) diffusion propagator 
and (cr(ri)<7(r 2 )) is the correlation function of the spins 
of LSDs. We can now estimate t w as 



1 



F 2 nS 2 



F 2 



1 



-7S 2 (l+fb) 2 



D t g 2 {l + F ) 



(67) 



where again the dominant contribution comes from LSDs 
of smallest size [n is the concentration of LSDs, see 
Eq. (fill)] . This should be compared with-the contribution 
of the usual Gaussian spin fluctuations^ 



2F 2 



T 



t* (1 + F )(2 + F ) g 



ln 5 (l + F ). 



(68) 



Again, up to the numerical factors, the contribution of 
LSDs dominates at temperatures lower than T* [given by 
Eq. (§)]. 

The actual crossover temperatures for different physi- 
cal quantities might differ by a factor of order In g(l+F ), 
but such difference is beyond the accuracy of our treat- 
ment. However, the discussed contribution of LSDs to 



VII. CONCLUSIONS 



In this paper we have considered the effect of disorder 
on magnetic properties of the ground state of a metal 
close to the Stoner instability. We have shown that even 
though on the mean field level the ground state of the 
metal is paramagnetic (1 + Fq > 0), there is non zero (ex- 
ponentially small) probability to form local spin droplets, 
i.e. domains of non zero spin polarization. The proba- 
bility to form a LSD is independent of its size R, thus 
LSDs of any size can appear. The total spin of the LSD 
is also independent of its size and obeys the distribution 
Eq. (|^), with the typical value S ~ (1 + i*b) _1 , which 
is large 1 <C S <C g, so that LSDs can be considered as 
classical spins. 

Considered as independent moments, LSDs contribute 
to the observables, changing the temperature dependence 
of both the paramagnetic susceptibility and the dephas- 
ing time at temperatures lower than certain cross-over 
temperature Eq. (^35|) . When T < T* , the dephasing time 
saturates to the temperature-independent value Eq. (|67|), 
while the susceptibility acquires the Curie-like 1/T de- 
pendence. 

Both the Curie susceptibility Eq. (|64|) and the dephas- 
ing time Eq. (^) were obtained in approximation of non- 
interacting LSDs. This approach is valid at temperatures 
larger than the typical value (^2|) of the interaction be- 
tween LSDs. Since the cross-over temperature T* is much 
larger than the interaction (|62|), there is a parametri- 
cally wide temperature regime [by our large parameter 
g(l + F ), see Eq. @], U t < T < T*, where the Curie 
behavior of the susceptibility and the saturation of the 
dephasing time can be observed. At smaller tempera- 
tures T < Ut, however, the LSDs can not be considered 
as non-interacting moments and the behavior of the sys- 
tem changes. Interaction of LSDs with each other or with 
itinerant electrons should lead either to screening of the 
local spins or to forming some spin glass state (due to 
the random sign of the interaction). Such regime was 
not considered in this paper.. 

We acknowledge helpful discussions with L.B. Ioffe and 
B. Spivak. We also acknowledge the hospitality of Theo- 
retische Physik III, Ruhr-Universitat Bochum, where this 
work was started. LA. is Packard research fellow. A.L. 
is supported by the NSF under Grant 9812340. 
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APPENDIX: 

Here we discuss the correlation functions that describe 
the mesoscopic fluctuations in the system, i.e. give the 
weight in the distribution Eq. (^Tj) of the random quan- 
tities F\ and B. Both quantities are the coefficients in 
the expansion of the thermodynamic potential Eq. (^) in 
powers of the spin density. Therefore the diagrams for 
Fi and B can be obtained by differentiating the diagram 
for the exact thermodynamic potential. To calculate the 
correlators one has to multiply the random quantities and 
then average over the disorder. As a result we get three 
different correlators, depicted diagrammatically on Figs. 
2-4. 

In terms of exact electronic Green's functions the ther- 
modynamic potential can be written as 



G R (e) 



(Al) 



where the Green's function is defined in the fluctuating 
field a, i.e. 



G R ( A \e) = \e-H-Ta±iO 



(A2) 



Here H is the Hamiltonian of the system, t j are the Pauli 
matrices. The symbol Tr denotes the trace over spin in- 
dices and the integration over all spatial cooridinates. 
For brevity we do not explicitly indicate that G R ( A \e) 
depend on spatial coordinates (due to disorder in H). 

The expansion (|^) can be achieved by taking varia- 
tional derivatives of the thermodynamic potential (Al) 
with respect to a. The second variational derivative of 
the thermodynamicpotential corresponds to the second 
order term in Eq. (0), 



s 2 n 



5a a (x 1 )da> 3 {x 2 ) 



de 



2ni 



xTr' G R T a G R T 



G A T a G A T 



(A3) 



Here the prime in Tr' means that there is no integration 
over the coordinates in the left-hand side (in this case x\ 
and x 2 ). The trace over the Pauli matrices gives (in the 
absence of spin-orbit coupling) Tr [t q t^] = 25 a/3 . Upon 
subtracting the average, the derivative ( |A3| ) is propor- 
tional to the non-local quantity i*i(:Ei., (x 2 ): 



Fi (x 1 ,x 2 ) 



5 2 n 



Sa a (x 1 )Saf 3 {x 2 ) 



S 2 fl 



x 5a a (x 1 )5al 3 (x2)' 



(A4) 



Similarly, the coefficient B in Eq. (|J) is given by the 
fourth varia tional derivative of the thermodynamic po- 
tential ( [AlD (only the average is zero in this case): 



G R T a G R T G R T >i G R T u 



_ G A r a. G A T G A rl XQA T u 



(A5) 



The trace over the Pauli matrices has two parts, 
Tr[r a T /3 T' i T !/ ] = 2[S a(i S^ -(S^S^ -S^S^t 1 )). The first 
part corresponds to the fourth order term in Eq. (|^), 
while the second part is the cross-product term, which 
disappears since we consider only the singlet fluctuations 
a =(0,0, a). 

We are now ready to calculate correlation functions of 
the mesoscopic fluctuations Fi and B. Averaging over 
disorder is performed in a standard manner (see Ref. @ 
for details). The correlator of the fluctuation Fi(xi,x 2 ) 

K FF [{x 3 },{yj}\ = 

K FF (x 1 ,x 2 ,y 1 ,y 2 ) = (F 1 (x 1 ,x 2 )F 1 (y 1 ,y 2 )) (A6) 

is given by a sum of four diagrams, one of which is de- 
picted on Fig. 2. 





FIG. 2. Diagrams for the averaged correlation function of 
the non-linear part of the paramagnetic susceptibility. 

All four diagrams contribute equally to the integral 
Eq. ([55]) (as will be clear from the explicit form of K FF ), 
therefore we shall proceed evaluating the contribution of 
the diagrams on Fig. 2. After averaging, each diagram 
corresponds to a product of four diffusion propagators 
(or diffusons) and four vertex blocks. 

The diffuson in momentum representation is given by 



V{u;q) 



1 



Dq 2 



where D is the diffusion constant, and q is a 2D mo- 
mentum vector. Each vertex block (see Fig. 4) contains 
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precisely the same factors 2-kvt 2 as each diffuson, thus 
in momentum representation the diagrams on Fig. 2 are 
expressed in terms of the integral 

o o 

K FF [{k j }] = J ^ J ^ Jd 2 QV(Lu;Q) 



To estimate the size of the fluctuation regions with non 
zero spin (LSD) we need the explicit dependence of the 
correlator Eq. (All) on the parameters of the problem. 
To do that we introduce a length scale R whi ch charac- 
terizes the size if the LSD and write Eq. (All) as 



xV(uj; Q - fei) V(oj; Q - fc 3 ) V{uj; Q-k x - k 4 ) 



x 5(ki + k 4 — k 2 — fc 3 ), 



(A7) 



where u> — e\ — e 2 . 

For the purposes of this paper it is more convenient to 
write the correlator Eq. (A6) in the coordinate represen- 
tation as (we have also evaluated one of the frequency 
integrals) 

oo 

/d 



xT>(lj; xi - y 2 ) V{oj; x 2 - yi) V(u; y 2 - x 2 ), (A8) 
where Xi and yi are 2D coordinate vectors and 
d 2 q 



V(u>;x) 



(2tt) : 



(A9) 



is the Fourier transform of the diffusion denominator to 
the position space. 

The remaining frequency integral can be evaluated by 
using the following integral representation for T>(ui;x). 
First, we represent the diffusion denominator as an inte- 
gral over an auxiliary variable 



— iuj + 



oo 

^— 2 = J dtexp[-t{-iuj + Dq 2 )]. 



The momentum integral in Eq. (A9) becomes Gaussian 
and we obtain 



1 



AttD J t 
o 



dt 



exp 



ILOt - 



UD 



(A10) 



Substituting the integral representation Eq. (A10) into 
the correlator Eq. (AS), we find the final expression for 
the cont ribu tion of the diagram on Fig. 2 to the correla- 
tor Eq. (|A6|) 



K FF[{Xj},{yj}} = 



1 



1 



dtidt 2 dt 3 dt4 

(2nvD) 2 4tt 4 J ht 2 t 3 t 4 
o 



exp 



ti 



t 2 



ta 



U 



(h + t 2 + 1 3 + t 4 ) 2 



(AH) 



K FF[{Xj},{yj}} 



9 2 R 4 



K 



FF 



(A12) 



where Kp p is the dimensionless counterpart of the inte- 
gral Eq. ([All]). 




FIG. 3. Typical diagram for the averaged correlation func- 
tion of (FiB). The rest of the diagrams are obtained by shift- 
ing positions of the spin vertices on the inside line relative to 
those on the outside line and interchanging coordinate indices 
similar to the diagrams on Fig. 2. There are 3 topologically 
different diagrams. 




FIG. 4. Typical diagram for the averaged correlation func- 
tion of (BB) . The rest of the diagrams are obtained by shift- 
ing positions of the spin vertices on the inside line relative to 
those on the outside line and interchanging coordinate indices 
similar to the diagrams on Fig. 2. There are 8 tpologically 
different diagrams. 
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The remaining correlator s in E q. ( |1 1| ) are constructed 
in the same manner as Eq. ( All ), the only difference be- 
ing the number of diffusons. The corresponding diagrams 
are given in Fig. 3 and 4. Again, the total number of 
diagrams is large, therefore we give expressions for the 
typical contributions depicted on Figs. 3 and 4. The rest 
of the diagrams are obtained by interchanging coordinate 
indices. 

The diagram on Fig. 3 contains six diffusons and six 
vertex blocks. Since each of the vertex blocks carries a 
factor of i, the overall sign of the di agram is negative. 
Using the integral representation Eq. (A10) for diffusion 
denominators, we hnd 



K F B[{xj},{Vj}] 



1 



1 



2(2vr) 4 (2-kD) 4 



n 

k=l 



dtk expj-7>a] 

k (E U) 2 

i=l 



(A13) 



T, 



FB 



(xi - y-i) 2 , (xi - va) 2 (x 2 - V2f 



tl ti tz 

(x 2 - yi) 2 . (yi - yi) 2 . (2/3 - 2/2) 2 



+ 



U 



ts 



The correlator on Fig. 4 can be written in the same 
way. It has positive sign, since it contains eight vertex 
blocs. 



Kbb [{xj},{yj}] 



1 



8(2tt) 4 (2ttD) 



n 

fc=l 



dt k exp[—T BB ] 



tk 



(A14) 



Tbb = 



(xi - y-i) 2 . (xi - yi) 2 (x 2 - ?/2) 5 



ti t<i t-$ 

jx2 -yi? {x 3 -y?,) 2 {x3 - y2? 

ti t$ tg 

(x 4 - yi) 2 (x 4 - yi) 2 

Altogether, these correlators can be combined in the 
matrix Eq. ( |ll| ) 

a _ ( K FF [{x 3 },{yj}} K F B[{xj},{yj}]\ , A1 - 
{K FB [{x 3 },{ yj }} K BB [{ Xj },{ yj }})- [ALb) 



Similarly to Eq. (A12), the dimensional analysis gives the 
size dependence in terms of the scale R 



ic = 



1 



9 2 R 4 



jpK FB jpK BB 




(A16) 



The sign of the correlator follows from the diagrams, so 
that the dimensionless integrals K are positive numbers. 
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